Essential state of the electromagnetic field and the double-slit experiment 
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A new class of generalized solutions related to the essential spectrum of linear Maxwell's equations 
is presented. The essential modes are given in terms of normalized singular Weyl's sequences, whose 
square represents Dirac's delta functions in spatial and angular frequency domains. The action 
integral associated with essential modes is well-defined. We claim that these modes represent the 
collapsed state of the electromagnetic field and, with some additional assumptions on the conserva- 
tion of action, are suitable for describing the double-slit experiment in accordance with the orthodox 
point of view. 



Wave motion plays a fundamental role in both clas- 
sical and quantum physics. Mathematically, the propa- 
gation of classical (physical) waves as well as evolution 
of quantum-mechanical wavefunctions (probability am- 
plitudes) are described by some kind of a linear wave 
equation. In quantum physics, however, wave motion is 
not everything. During the measurement wavefunction 
undergoes what is known as wavefunction collapse. This 
collapse is one of the phenomenological postulates of the 
orthodox quantum theory. It remains a postulate, since 
such collapses are believed to be mathematically incom- 
patible with the linear wave motion in general. Proba- 
bly, the most well-known and at the same time the most 
paradoxical case is the double-slit experiment with pho- 
tons, where a spatially extended wave collapses into what 
seems to be a truly random point at the detector screen. 
A large collection of such points, however, accumulates 
into a classical interference pattern. 

Although, the standard quantum optics operates 
within the postulates of the orthodox theory and pro- 
vides excellent predictions about everything that has to 
do with photons, many find the pragmatic orthodox view- 
point difficult to accept and look for some kind of a pro- 
cess, which would resemble the collapse. There are two 
possibilities in this respect. One is to find the explanation 
for collapse within the quantum physics itself. For exam- 
ple, the decoherence approach [jj] considers the interfer- 
ence of partially coherent wavefunctions entangled with 
the environment. Although this provides a better under- 
standing of the quantum-classical transition, supporters 
of the decoherence approach tend to dismiss the ortho- 
dox collapse as an awkward and unnecessary theoretical 
construct. The other possibility is to find the classical 
solution of a wave equation which resembles the collapse, 
as done in the "dynamical reduction" program 2]. This 
program is largely based on the ideas from classical non- 
linear wave motion. However, nonlinear focusing is never 
perfect, so that the idealistic orthodox collapse is dis- 



missed within this approach as well. Moreover, nonlinear 
wave propagation is a causal and deterministic process. 
Therefore, in addition to an artificial nonlinearity in the 
modified Schrodinger equation the dynamical reduction 
program introduces a stochastic term as well. 

Here, somewhat contrary to the growing consensus, we 
would like to defend the orthodox point of view, namely, 
the idea of collapse. Ontological (interpretational) prob- 
lems aside, the main shortcoming of this idea is its pos- 
tulated nature. Of course, we do not pretend here to 
know the general solution of this problem on an abstract 
level. However, in the case of the electromagnetic wave 
motion, we do seem to have found a unified description 
of the propagating and collapsed electromagnetic field, 
so that both causal electromagnetic waves and the ran- 
dom points at the detector come out as solutions of linear 
Maxwell's equations. This is the subject of our Letter. It 
is difficult to decide whether our approach is classical or 
quantum-mechanical. Probably, it is neither. And, may 
be, it is only natural, since the problem lies somewhere 
at the vague border between quantum and classical. 

To accommodate for collapse one needs a perfectly lo- 
calized natural mode of the electromagnetic field, which 
solves linear Maxwell's equations in a certain, possibly 
generalized, sense. Elsewhere we show that the required 
mode is related to the essential spatial spectrum of the 
electromagnetic scattering operator |3|, and is given in 
terms of a singular Weyl's sequence jj]. The action in- 
tegral and other quadratic physical quantities involving 
the electromagnetic field are well-defined for this spa- 
tially localized mode. Thus, despite its generalized na- 
ture, we could not dismiss the obtained singular mode as 
a purely mathematical artifact and looked for its physical 
interpretation. It turns out that our work is closely re- 
lated to experiments in photonics and plasmonics, which 
both exploit unusual resonant states of the electromag- 
netic field. Recently we have uncovered a link between 
the plasmonic resonance and the essential spectrum of 



2 



the electromagnetic scattering operator 

However, a spatially localized mode, even a perfectly 
localized one, is not quite enough. Although, the con- 
nection is not immediately obvious, states of this kind 
are already embedded in the modern mathematical ap- 
proach to quantum mechanics known under the collective 
name of Rigged Hilbert Space (RHS) |(|. Apparently, 
the generalized eigenf unctions of the essential spectrum 
mentioned above do not belong to the usual complete 
Hilbert space Q. The RHS approach allows, in princi- 
ple, the quantum evolution to escape beyond the Hilbert 
space. This idea, however, is still at odds with another 
unpopular part of the orthodox interpretation. Namely, 
suggestion that collapse happens instantaneously and is 
not the result of a causal evolution. Again, the general 
consensus is to replace the instantaneous orthodox quan- 
tum "jumps" with a possibly very fast, but continuous 
process, e.g. evolution of a quantum system into a deco- 
herent state. 

Retaining the idea of a spatially localized resonant 
state, we suggest the following paradigm shift. First 
of all, we know that the discrete spatial spectrum of 
the Maxwell operator corresponds to delocalized (global) 
eigenmodes of the electromagnetic field, either con- 
fined within a resonator our rapidly decaying at infinity. 
Whereas, as we show in and here, the essential spa- 
tial spectrum corresponds to perfectly localized and nor- 
malizable generalized modes. By analogy we introduce 
the discrete temporal spectrum (spectrum of the time- 
derivative operator) corresponding to finite-duration pro- 
cesses, which have both the beginning and the end, e.g. 
a photon is first radiated and then absorbed. The essen- 
tial temporal spectrum corresponds to generalized quasi- 
harmonic "processes" , properly normalized and perfectly 
localized in the angular frequency domain. Combining 
these observations we propose to distinguish between the 
"normal" and the "essential" spatio-temporal states of the 
field. 

To be more precise mathematically, in the electromag- 
netic case we are talking about the "Weyl" rather than 
essential spectrum, as we apply the Weyl criterion Q to 
the differential Maxwell's equations, and the operator in 
question is not selfadjoint. This however is not impor- 
tant for the following two reasons. First of all, the spatial 
essential (singular) modes we encounter here are, in fact, 
proven to be the generalized modes of the essential spec- 
trum of the equivalent non-selfadjoint and non-normal 
integral operator 0,0]. Secondly, for our purposes the 
generalized (and normalized) nature of the modes cap- 
tured by the Weyl criterion is much more important than 
the actual name of the spectrum they belong to. Con- 
sider the Maxwell equations 

-V x H + d t D = -J, 
V x E + d t B = 0, ^ ' 

with linear isotropic constitutive relations D(x, t) = 



e(x, t) * E(x, t), B(x, t) = t) * H(x, t), where e(x, t) 
and t) are some continuous functions with finite spa- 
tial support, and star denotes the temporal convolution. 
In the matrix form these equations can be written as 



(2) 



In operator notation we simply write MX = Y, where M 
is the Maxwell operator, and X and Y are six-vectors. 
We introduce the Hilbert space with the following norm: 
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According to the Weyl's definition of the spectrum of a 
selfadjoint operator M one needs to find the sequence of 
functions {F„} with the following properties: 



= 1, 



lim 

n — >oo 



MF„-AFJ| =0 



(4) 



If the Weyl sequence that satisfies these equations for 
some A has no convergent subsequence on the pertaining 
Hilbert space, then such A is in the essential spectrum. 
If for some other A the Weyl sequence does converge to a 
function from the Hilbert space, then this function is an 
eigenfunction, and A is an eigenvalue of M. With some 
reservations, which we mention earlier, we shall apply 
this definition to the Maxwell operator. The singular 
modes of the volume integral operator found in are 
the following vector-valued functions: 
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where x, x c G R 3 and a > is the sequence parame- 
ter. The properties of these functions are summarized in 
Theorem 2.1 of The most important are the following 
two: 

f |*(a,x,x c )| 2 dx=l, (6) 
lim f /(x)|*(a,x,x c )| 2 dx = /(x c ), (7) 

which is why we refer to singular modes as the square 
root of the delta function. As we intend to take into ac- 
count the temporal spectrum, we shall consider a similar 
sequence of scalar functions: 

<j)(/3,t,uj c ) = 

iA/27T- 1 /4 yS -3/4 texp (^-^t 2 - lLU c ?j . 

The Fourier transform of these functions is 

^vr-^/V-c^exp (-i\co - Uc \ 2 ^ , 



(8) 



(9) 
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and, as can be verified using the integration techniques 
developed in Q, it is normalized and also represents the 
square root of the Dirac delta function, but now in lu- 
domain, i.e., 



<p{f3,0J,uj c 



dLU = 1, 



(10) 





F m = 










0* 



f°° „ 2 

lim / g(uj) 4>((3,u,u c ) du; = g(u c ). (11) 

The product of functions (0 and © can be used to gen- 
erate the essential modes of the complete Maxwell's op- 
erator, i.e., 



(12) 



Indeed, substituting F c and F m in (0} we eventually ar- 
rive at 

lim ||MF°- A°F C || 2 = |A c + z^ c e(x c ,^ c )| 2 , (13) 
Km ||MF m - A m F m || 2 = |A m + iuj c fi(x c , u; c )\ 2 , (14) 

a,/3 — >oo 

meaning that — iw c e(x c , lo c ) and — ic<j c /}(x c , u> c ), where 
u> c 6 R and x c G R 3 , are in the essential (or Weyl's) 
spectrum of the Maxwell operator. 

Mathematical suitability of the essential state repre- 
sented by one or more essential modes for describing the 
process of photon detection should be apparent by now. 
Indeed, the photon is detected at a single point of the 
detector-screen, which is x c in our case, but also trans- 
fers a discrete portion of energy, which in our case may be 
associated with the angular frequency lo c . This, however, 
is a tricky question, as we do not know, if any energy can 
be assigned to the essential mode at all. 

What are the mathematical and physical conditions 
on the excitation of essential modes? From the Maxwell 
equations we deduce that F c and F m can only be excited 
by either induced or external currents with generalized 
components of the form — Sod t F c or —fi dtF m . This is 
applicable to both macroscopic and microscopic electro- 
dynamics. In the macroscopic linear case with induced 
currents J ind = ft[(e*)-e ]E and K ind = dt[(p*) ~Mo]H 
this leads to the condition on the essential resonance, i.e. 
A c = — iw c e(x c , ui c ) — and A m = —iui c jj,(x c ,uj c ) = 0. 

We notice another important similarity with the ortho- 
dox collapse, as essential states cannot causally evolve 
out of an initially normal state, even in principle. If 
this could happen, then there would exist a sequence of 
Hilbert space currents causally related to the sequence 
of fields, both converging to an essential state. However, 
this is not possible by the very definition of the essential 
state. Hence, the normal — > essential state transforma- 
tion is a non-causal "jump-like" event postulated by the 
orthodox theory. 



Thus we cannot say much about the normal — > essen- 
tial state transformation at the level of the fields. If, 
however, the field is squared and integrated, then both 
the normal and the essential states can be treated on 
equal footing. Note that a higher power (e.g. cube) of 
an essential mode gives an infinity, while the lower power 
produces a zero 8]. We shall approach the problem of 
the transformation in the following way. Suppose that 
the normal state does not evolve into, but is simply re- 
placed by an essential state. An equation describing such 
an event can take various forms and as we show below 
is equivalent to the conservation of the electromagnetic- 
field action in the presence of essential modes. 

Assuming that the normal — » essential replacement 
transformation takes place, the practical problems one 
should try to solve in relation to the double-slit ex- 
periment are: the spatial locations of the essential 
modes, their angular frequencies, and the precise amount 
of action they carry (action is the most elementary 
squared-and- integrated quantity involving fields). But 
first we note an important degree of freedom, a kind 
of non-uniqueness, associated with the essential modes. 
Namely, condition for the essential resonance A c = 
— iw c e(x c ,w c ) = 0, if satisfied at all, will, in general, be 
satisfied not just at a single point in space, but along 
a hyper-surface or even a hyper-volume. If a zero per- 
mittivity sounds like something suspiciously unphysical, 
one is advised to think of essential induced currents in- 
stead, without any reference to the macroscopic consti- 
tutive relations. In any case, we shall assume that the 
measurement apparatus is prepared in such a way that 
the conditions on the excitation of essential modes are 
satisfied at the entire detector screen and for some range 
of angular frequencies. From this it follows that certain 
expressions containing a single squared and integrated 
essential mode will be invariant with respect to the lo- 
cation x c and frequency lu c as long as this location is 
somewhere at the detector screen and the frequency is in 
the prescribed range. 

Indeed, let E nr (x, t) and E cs (x, t) denote the elec- 
tric fields in the normal and essential states, respec- 
tively, undergoing the replacement transformation. For 
£o | E os 1 2 = C|F C | 2 consider the following basic conserva- 
tion law: 



n 


E nr (x,w) 







C lim 

OL,0 — K> 



dx duj 



4>(0,w, w c )*(a,x,x c 



(15) 



dxdw = C, 



where C is a fixed dimensional constant, which can 
and, probably, should be chosen as H - the Plank con- 
stant, since it represents the elementary portion of ac- 
tion. To properly justify the above law one should, per- 
haps, look into the electromagnetic field part of the la- 
grangian /io|H| 2 — £o|E| 2 (which is the same in classi- 
cal and quantum electrodynamics). For instance, we see 
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that the simultaneous addition of £o|E os | 2 = C|F°| 2 and 
/io|H os | 2 = C|F m | 2 does not change the value of the ac- 
tion integral. Another possibility is to consider not only 
the smooth variations of the action integral as we do 
when we recover the classical equations of motion, in- 
cluding the Maxwell equations, but also its "essential" 
variations. This, however, is beyond the scope of the 
present Letter. 

Obviously, the basic conservation law l|15|) does not 
depend on x c and lo c . However, if we are to describe 
the double-slit experiment, the location and frequency of 
a number of modes has to be distributed in accordance 
with the intensity (squared amplitude) of the field in the 
normal state - the Born statistical postulate. We can 
show that the Born postulate is equivalent to a more 
general conservation law associated with the normal- 
essential state transformation. Let each essential mode 
F c carry a finite portion of action as in equation (JTSJ. 
Hence, position x c and frequency lu c of each particular 
essential mode are, in general, not defined (arbitrary). 
We could say that they are "random" , but this assump- 
tion does not seem to be necessary. Further, for any func- 
tion w(x, w), let the essential states satisfy the following 
generalized conservation law: 



So 



so 
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it(x, U>) 


E nr (x,w) 


I -oo Jxei 


I 3 
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w(x, to) 


E os (x, uj) 


/ -oo Jx£l 


i 3 





dxdw = 



dx duj, 



(16) 



where expression on the right contains all the necessary 
limits over Weyl's sequences. A similar law must hold 
for the magnetic fields. Our goal is to determine the dis- 
tributions of x c and uj c for the resulting essential state. 
Suppose that the normal state is transformed into N es- 
sential modes, then 



sq 
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7i(x, Uj) 


1 -oo JxGl 


I 3 



dxdcj = 



lim NC 

a,{3 — >oo 



oo «/x(E 



1i(x, Uj) 

0(/3,w,o; c )*(q!,x,x c ) 

it(x, w)/(x)g(w) dx duj, 



oo J x c G 



dx r duj r dx duj 



NC 



f(x c )g(uj c )x 



(17) 



where /(x c )dx c and g(uj c )duj c are the fractions of essen- 
tial states with parameters within the corresponding sub- 
volume and interval, i.e., /(x c ) and g{ui c ) are the distri- 
butions to be found. For arbitrary functions m(x, uj) this 
equality holds if and only if 



so 
NC 



E nr (x,w) 



(18) 



In other words, although for each particular essential 
mode parameters x c and lu c are not defined, their dis- 
tribution is equal to the normalized squared amplitude 



of the electromagnetic field in the normal state, i.e. for 
large N the essential modes should form an interference 
pattern. 



The electromagnetic field can, apparently, exist in two 
distinct states. One, related to the discrete space-time 
spectrum, is the causal state, i.e., the classical electro- 
magnetic wave. The other state, related to the essen- 
tial space-time spectrum, represents the collapse of the 
electromagnetic field into generalized electric/magnetic 
modes with undefined position and angular frequency. 
These modes are given in terms of singular Weyl se- 
quences that do not converge (in norm) to any function 
on the Hilbert space. The two states cannot causally 
evolve into each other. However, they may, under right 
circumstances (e.g. essential resonance), undergo a re- 
placement transformation obeying certain conservation 
laws. In particular, we have made two plausible assump- 
tions about the normal — > essential state transformation. 
First, each essential mode can carry (absorb) only a fi- 
nite portion of action, which, in fact, is an extension 
of the Plank law. Second, these transformations must 
obey a generalized conservation law. Under these as- 
sumptions the Born statistical postulate is uniquely re- 
covered. Thus, it would be interesting to find a proper 
justification for the mentioned generalized conservation 
law. If we are right about the nature of collapse, then 
there should be a higher-level equation which governs the 
rate of collapse. Such an equation would, probably, in- 
volve action integrals and describe the overall ratio of 
normal and essential states and the change of this ratio 
with time. 



[8] 



Electronic address: n.budko@ewi.tudelft.nl 
This research is supported by NWO (The Netherlands) 
W. H. Zurek, Reviews of Modern Physics, 75, 715 (2003). 
A. Bassi and G. C. Ghirardi, Physics Reports, 379, No. 5, 
257-426, (2003). 

N. V. Budko and A. B. Samokhin, SIAM J. Sci. Comput., 
28, No. 2, pp. 682-700, (2006). 

N. V. Budko and A. B. Samokhin, (under review) SIAM J. 
App. Math., also to be found at arXiv, physics/0606137 
(2006). 

N. V. Budko and A. B. Samokhin, Phys. Rev. Lett., 96, 
023904, (2006). 

A. Bohm, Phys. Rev. A 60, 861876 (1999). 
P. D. Hislop and I. M. Sigal, Introduction to Spectral 
Theory: With Applications to Schrddinger Operators, 
(Springer- Verlag, New York, 1996). 

J. F. Colombeau, Multiplication of Distributions: A Tool 
In Mathematics, Numerical Engineering and Theoretical 
Physics, (Springer- Verlag, Berlin, 1992). 



